Stability of compressible reduced magnetohydrodynamic equilibria 
— analogy with magnetorotational instability 
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Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) 
model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared 
with those obtained for ideal magnetohydrodynamics (MHD) from the classical SW criterion. An 
identification of the terms in the second variation of the free energy functional for CRMHD with 
those of SW is made: two destabilizing effects present for CRMHD turn out to correspond to the 
kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending 
and compressibility are also identified in the reduced model. Also, using the EC method, stability 
conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a 
reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. 
Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the 
condition for magnetorotational instability (MRI) for the latter model, corresponds to the condition 
for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of 
the curvature term and pressure gradient, respectively. Using the EC method, stability conditions 
for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, 
are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to 
the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, 
the Rayleigh hydrodynamical stability condition is recovered. 



I. INTRODUCTION 

The investigation of instabilities is of paramount im- 
portance for fusion and astrophysical plasmas. For fusion 
plasmas, detecting the main instabilities is necessary to 
learn how plasma characteristics such as current, den- 
sity, and temperature should be optimized in order to 
improve plasma confinement. For astrophysical plasmas, 
the identification of instability mechanisms is essential 
for understanding phenomena such as solar flares, coro- 
nal mass ejections, transport of angular momentum in 
accretion disks, and several other major observable phe- 
nomena. In the framework of the fluid description of 
plasmas, the use of the energy principle is one of the pos- 
sible ways to investigate linear stability in the limit where 
dissipative processes are negligible. 

In their early work, Bernstein et al. pQ provided the 
energy principle for the investigation of stability in ideal 
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MHD. According to such principle, an ideal MHD equi- 
librium is stable if and only if the quadratic form SW 
(more aptly called S 2 W, but we follow tradition), rep- 
resenting the second variation of the potential energy, 
induced by variations of the field variables, is positive or 
zero, for all the allowable Lagrangian displacements of 
the fluid. Further extensions of the energy principle, and 
application of it to particular geometries have been de- 
veloped in the subsequent years (see, e.g., Ref. 0), and 
methods based on the energy principle are now standard 
for investigating plasma stability. 

Alternatively, the EC method, an Eulerian stability 
method with an early antecedent in plasma physics [3], 
evolved into a formalized and systematic procedure be- 
cause of the discovery of the noncanonical (Eulerian) 
Hamiltonian description of ideal MHD in Ref. [I]. Of 
particular relevance to the present paper are a large num- 
ber of works that followed [4] detailing the noncanonical 
Hamiltonian formulation of reduced models, viz., reduced 
MHD [31 [5] , the Charney-Hasegawa-Mima equation [7] , 
the four-field model for tokamak dynamics [8], models 
for collisionless reconnection [9l [10] , hybrid fluid-kinetic 
models and models for many other systems. (See 
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[HHH for review.) 

Hamiltonian structure provides a natural means of ob- 
taining sufficient conditions for stability, and this is gen- 
eralized for the Hamiltonian structure of fluid models 
for plasmas that are defined on an infinite-dimensional 
phase-space and are noncanonical. The latter property 
gives rise to the existence of a special class of invariants 
for the model: the Casimirs. These are functionals of the 
observables, that commute, through the Poisson bracket 
of the system, with any other functional of the observ- 
ables. The identification of the Casimirs for a system, 
provide the basis for the EC method (see, e.g., O fToT - 
IT7] for a review). 

According to the EC method, stability is attained if the 
second variation 8 2 F, where F is a free energy functional 
obtained as a linear combination of the Hamiltonian and 
the Casimirs, has a definite sign, when evaluated at the 
equilibrium of interest. This method can be applied to 
the huge class of plasma models possessing noncanonical 
Hamiltonian structure. It provides sufficient conditions 
for energy stability, a strong form of stability that implies 
linear stability (while the converse is in general not true) 
and almost nonlinear stability. 

Given the above considerations, the question naturally 
arises of whether there is a connection between the stabil- 
ity results obtained from the EC method and those of the 
SW energy principle. One goal of this paper is to explic- 
itly trace a connection between the two methods in the 
context of two reduced fluid models: the CRMHD model 
and the accretion disk model introduced in Ref. |18) . in 
their dissipationless limits. 

CRMHD [19 , was conceived to investigate tokamak 
dynamics in a simplified geometry, using a large aspect- 
ratio ordering. It generalizes reduced MHD by taking 
into account finite compressibility, which involves a cou- 
pling between parallel flow and pressure evolution. It 
accounts also for effects due to magnetic field curvature. 
Its Hamiltonian structure was derived in Ref. [8]. 

The accretion disk model, on the other hand, describes 
locally an incompressible rotating plasma in the presence 
of an imposed azimuthal sheared flow and a poloidal field. 
It has been derived, via an asymptotic expansion, in order 
to investigate, in a simplified setting, the MRI [2"Ull23| . 
which is believed to play a key role in the transport of 
angular momentum in accretion disks [21] ■ The Hamil- 
tonian structure for this model, obtained here, was not 
previously known. 

In this paper, we carry out an analysis in the case of 
CRMHD, first. For a wide class of equilibria, we are able 
to trace the connection between the two methods. In 
particular, the terms that can lead to indefiniteness in 
the sign of S 2 F, actually correspond to the terms associ- 
ated with the kink and interchange instabilities in SW. 
We identify also a correspondence between the stabiliz- 
ing terms that appear in the expressions obtained from 
the two methods. 

The analysis in the case of the accretion disk model fol- 
lows essentially as a by-product of an identification rule 



that permits us to map one model into the other. This 
provides us with a convenient framework for achieving 
another objective of this article, which is the identifica- 
tion of analogies between instabilities relevant for toka- 
maks and those for accretion disks. By taking advantage 
of the mapping, we are able to show that the MRI model 
possesses a Hamiltonian structure. We apply then the EC 
method to this model and derive energy stability condi- 
tions. In addition to recovering the expected condition 
for MRI, we are able to see that such instability in the 
accretion disk model, plays the role of the interchange 
instability in CRMHD. Thanks to the unifying frame- 
work provided by the EC method and to the connection 
with the energy principle, we can establish further analo- 
gies, from which emerge quite strong similarities between 
the two physical systems, providing an example of uni- 
versality between instabilities in fusion and astrophysical 
plasmas. 

The article is organized as follows. In Sec. [TTJ we re- 
call the CRMHD model, briefly review the noncanoni- 
cal Hamiltonian description for continuous media, and 
then derive stability conditions for CRMHD using the 
EC method. Next, stabilizing and destabilizing terms of 
the EC method are compared to those of the SW energy 
principle by mapping Lagrangian variations to Eulerian. 
In Sec. |III| we present the formal correspondence between 
CRMHD and the MRI model. Section [TV] is devoted to 
the stability analysis of the MRI model and to the iden- 
tification of the similarities with CRMHD. Finally, we 
conclude in Sec. |V] 

II. CRMHD 

In a Cartesian coordinate system (x, y, z) as depicted 
in Fig. [T] the equations for CRMHD are given by 

^ + [ < i ) a)] !/ + [VV,4-2b>% = 0. (1) 

^ + [<t>,i>] y = 0, (2) 
dv 

— + [M v + Mv = o, (3) 

^ + [<My + /%,V>L-2/3[Mk/ = 0. (4) 

Equations ([lJ-Q , describe the dynamics of a plasma in 
a magnetic field B = eVipx z+^l + ex)^ 1 + eb]z + 0(e 2 ), 
with e = a/Ro <C 1 being the ratio of the characteristic 
length scale of the poloidal plane, a, with the major ra- 
dius of the magnetic axis, Rq, and b accounting for dia- 
magnetic corrections. The fields v and p represent the 
parallel ion velocity and the electron pressure, respec- 
tively, whereas ui — V 2 </> with <fi being the electrostatic 
potential. All the fields of the model are independent 
of the coordinate z. The quantity h = x accounts for 
the effects of magnetic field curvature, while the param- 
eter /3 is the ratio of the background electron fluid pres- 
sure and the magnetic pressure. Finally, the quantity 
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[/, g] y = z ■ V/ x Vg is the canonical bracket, with re- 
spect to the variables x and y, between two functions / 
and g. 





FIG. 1: Sketches of the geometry and local coordinate sys- 
tem for the accretion disk model (left) and for CRMHD for a 
tokamak (right). 

Equation is the equation for the evolution of the z 
component of the vorticity, ui, associated with the E x B 
flow, while Eq. ^ expresses the frozen-in condition for 
the poloidal magnetic flux tp. Equation ([3| governs the 
parallel ion dynamics and Eq. Q determines the electron 
pressure evolution. 



A. CRMHD Hamiltonian structure 

The CRMHD system of Eqs. Q-Q is known to pos- 
sess a noncanonical Hamiltonian structure [25] , which 
we review here. To this end we first recall a few basic 
notions about infinite-dimensional Hamiltonian systems 
(for further details see, e.g., Ref. [13]). 

Given a field theory xH x ) = F l {x 1 i'" iX N )> f° r N 
fields (% 1 (x), • • • , x (x)), evolving according to the vec- 
tor fields F l , we say that such system is a Hamiltonian 
system if a Hamiltonian functional -H^xS 1 ' ■ X ) and a 
Poisson bracket {, } exist, such that the field equations 
can be written in the form 



X 1 = {X\H}, 



i = l,---,N, 



(5) 



where the Poisson bracket must be bilinear, antisymmet- 
ric, and satisfy the Leibniz and Jacobi identities. 

For continuous media, such as plasmas described in 
terms of Eulerian variables, the case for CRMHD, Pois- 
son brackets typically possess a noncanonical structure. 
In these cases there exist Casimir functionals, that is, 
functionals C(x 1 , • • • , X N ) that satisfy {C, F} — 0, for all 
functionals F(\ 1 ,--- ,X )■ Because, C — {C,H} = 0, 
for any Hamiltonian, Casimirs are invariants of the dy- 
namics. 

For the CRMHD, the Hamiltonian functional is given 
by 



H 



dxdy(\W(f>\ 2 + v 2 + \Wtp\ 2 +p 2 /f3). (6) 



All the terms in the Hamiltonian are amenable to a phys- 
ical interpretation: the first two terms represent the ki- 
netic energy associated with the E x B flow and with the 



toroidal flow, respectively. The third term is the poloidal 
magnetic energy, whereas the last term accounts for the 
internal energy of the electron fluid. 

The noncanonical Poisson bracket for CRMHD, on the 
other hand, is given by 



{F, G} = fdxdy(u[F m G u ] y 



v([F u ,G v ] y + [F v ,G u ] y ) 

(p + 2ph)([F m G p ] y + [F p ,G u ] y ) 

i>([Fu,G^]y + [F^jG^y 
— f3[Fp, G V ]y — (3[F V , Gp\y) 



(7) 



where subscripts of the functionals indicate functional 
derivatives. The operation ([7| can be shown to satisfy 
all the properties of a Poisson bracket, provided that sur- 
face terms vanish when integrating by parts. This can be 
accomplished, for instance, on a bounded domain, im- 
posing that the fields and the functional derivatives are 
periodic in the y direction and satisfy vanishing Dirichlct 
boundary conditions along x. The Hamiltonian Q and 
the bracket Q, combined according to (pH), yield namely 
the CRMHD system Q-(4j). Casimir functionals of the 
bracket ([7|, on the other hand, are given by 



Ci = 


j dxdy J-(i{j), 




C-2 = 


j dxdy v Af(ip), 




c 3 = 


j dxdy£(Tp)(p/l3 


+ 2h), 


c 4 = 


j dxdy (uj Q{ip) — 


v Q'ty) (p//3 + 2ft)) , 



where J-, AC, C and Q are arbitrary functions, and the 
prime denotes derivative with respect to the argument 
of the function. This set of Casimirs is close to but 
not equivalent to that for unreduced MHD with various 



symmetries [55] HZ? 
invariants of Eqs 



Because of the infinite number of 
8j) the dynamics of CRMHD is con- 
strained. In particular, the existence of the family of 
Casimirs C\ reflects the constraint of the frozen-in con- 
dition for the poloidal magnetic flux. The families Ci and 
C4, on the other hand, are remnants of the cross- helicity 
/ d 3 i v ■ B, which is a Casimir for isentropic ideal MHD 
(see Ref. [25]). whereas C3 derives from the conservation 
of magnetic helicity of ideal MHD 25 . Such relation- 
ships with the Casimirs of full ideal MHD make use of 
the proportionality between p and b, which is imposed 
when deriving CRMHD from the fluid plasma descrip- 
tion (see, e.g., Ref. [2T)]). 



B. CRMHD equilibria and energy stability 

One of the advantages provided by the Hamiltonian 
structure of a model is the possibility of using specific, 
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well-developed methods for Hamiltonian systems; in par- 
ticular, methods that have been developed for determin- 
ing stability conditions. For noncanonical systems, such 
as CRMHD, the EC method is particularly convenient, 
because it can take advantage of the infinite number of 
Casimirs possessed by the model. The EC method ap- 
plies to equilibria that are extremals of a free energy func- 
tional F, which is a linear combination of the Hamilto- 
nian and the Casimirs. For such equilibria, the method 
can provide sufficient conditions for energy stability - 
stability that arises from the existence of a constant of 
motion for the system, whose second variation has a def- 
inite sign, when evaluated at the equilibrium of interest. 
Because, the Hamiltonian and the Casimirs are constants 
of motion, then clearly F is also a constant of motion. 
Therefore, for noncanonical Hamiltonian systems the free 
energy functional provides a natural constant of motion 
from which one can determine energy stability. We re- 
mark that this type of stability implies linear stability, 
but the converse is not true in general. Also, we remark 
that this kind of stability almost assures nonlinear sta- 
bility, but in general nonlinear stability does not imply 
linear stability. 

The EC method evolved out of early plasma work [3] 
and there is now considerable prior application to plasma 
physics problems, as can be found in Refs. [T5HT71 15UW32) . 

In order to apply the EC method to CRMHD, we define 
first the free energy functional, which for CRMHD, in its 
most general form, reads 



F = H + C\+C 2 + C 3 + C 4 . 



(9) 



Equilibria to which the EC method can be applied are 
extremals of ([9]), i.e., those at which the first variation 



8F 



dxdy 



V 2 V> + Fty) + v M'W) + {p/p + 2h) C'ty) 



+ ojG' O) - v Q"($) (p/P + 2h) J 5^ 
+ (y +Af(tp)-g'(i>)(p/p + 2h)yv 

+ (p + C{tp)-v Q'(S>j)5p/P 



(10) 



vanishes. Because of the arbitrariness of the variations 
Slu, Sip, Sp and Sv, such equilibria must be the solutions 
for <j>, ip, p and v, of the system 

<t> = (ii) 

Vfy = F'[ilj) + vAr'{ilj) + {p/P + 2h)£'{ip) 

+wG\ii)-vG"{il))(jp/P + 2h), (12) 
v = -Af(4>) + g'(i>)(p/f3 + 2h), (13) 
p = -£(V0+«S'(V>), (14) 
for choices of G, F, JV and C. 



We remark first that the equilibrium equation (111 ex- 



flow is a poloidal magnetic flux function. In the follow- 
ing, however, we will consider only the case G = 0, cor- 
responding to the absence of poloidal flow, which is a 
standard case treated in linear stability analyses of ideal 
MHD. However, with the following transformation to a 
dynaflux function X' 



di(;Jl-G' 2 W, 



(15) 



which is valid for sub-Alfvenic flows (G' 2 < 1), equilib- 
ria with poloidal flow can be mapped into ones without. 
This transformation was first noted in Ref. [33] and gen- 
eralized in Ref. [31]. We will not pursue this further 
here, but emphasize that for every equilibrium without 
poloidal flow there exist a family of equilibria determined 
by G with poloidal flow. 



If Q = 0, then Eqs. (13) and (14) imply that toroidal 



flow and pressure are constant on poloidal magnetic sur- 
faces. Inserting (13) and (14) into (12), yields 



V 2 V = F'{4,) - v\i>)v{>iP) - p'(V0(KV>)//3 + 2/i) , (16) 

which can be seen as a generalized Grad-Shafranov equa- 
tion for the magnetic flux function ip. Once this partial 
differential equation is solved for tp, the corresponding 
solutions for the fields p and v can be readily obtained 
from the specific choices of the free functions C and Af. 

When investigating stability for a reduced model, such 
as CRMHD, it is natural to ask how the stability proper- 
ties of the reduced model square with those of the com- 
plete MHD equations. This can help in interpreting the 
stability conditions for the reduced model and in seeing 
what features have possibly been lost because of the re- 
duction. In the following we compare the stability analy- 
sis for CRMHD carried out by means of the EC method, 
with the classical energy principle analysis jTj, based on 
the 5W criterion, for the complete ideal MHD system 
with no equilibrium flow. 

In order to facilitate this comparison, we consider first 
the special case where Af = 0. Together with G = 0, this 
assumption implies static equilibrium for CRMHD, i.e., 
equilibria with neither poloidal nor toroidal flow. For this 
case, the second variation of the free energy functional 
reads 



8 F 



dxdy \VS<t)\ 2 + \VSip\ 2 + \Sv\ 2 + \Sp\ 2 / '/3 

{F"(^) + (p/p+2h)c"mm 2 

2C'{ip)5p5^/P 



(17) 



presses the condition that, at equilibrium, the poloidal Using the relation p(ip) = —C(ip) and rearranging the 
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terms, one can rewrite S 2 F as follows: 



S 2 F 



dxdy 



\V8(p\ 2 + \Sv\ 2 + \V6ip\ 2 



Sp — 



x Bp • Vp 



H\ - 



x B, 

Bl 



2 (z x Bp ■ Vg) (z x Bp ■ V/t) a 



£ 2 



B 2 

v 



Vh |^|- 



(18) 



where B p = VV> x z denotes the poloidal equilibrium 
magnetic field and j|| = — V 2 ip the equilibrium parallel 
current density (apart from the rescaling factor e). 

According to the EC method, equilibria for which S 2 F 
has a definite sign, are energy stable. Analyzing the ex- 
pression (18), one sees that the first four terms are pos- 
itive definite, whereas the last two terms can be posi- 
tive or negative and, consequently, could give rise to in- 
stability. The first of these two terms depends on the 
current density gradient. In the absence of magnetic 
field curvature effects (h = 0), this term can be writ- 
ten as — J dxdy\5ip\ 2 dj\\/ dip. Therefore, a current profile 
that is monotonically decreasing in ip has a stabilizing ef- 
fect. This effect had already been identified by means of 
the EC method in the context of incompressible reduced 
MHD [Hj. For CRMHD, the parallel current density is 
no longer a magnetic flux function, because of the pres- 
ence of curvature, which will then play a role in determin- 
ing the stability condition. The second term of indefinite 
sign, on the other hand, is stabilizing if the gradients of 
the equilibrium pressure and of the magnetic curvature 
terms point in the same direction. 

We show in the following, that the two terms of indef- 



inite sign in ( 18 1 are related to the kink and interchange 



instabilities for ideal MHD. 

For comparison, consider the usual energy principle 
analysis for MHD. Recall, the expression for 6W in the 
form presented in Ref. [35] : 



8W = - I drill/ 



Q 



B 2 



B 



J 



^xB).Q 
(19) 



In (19 1 £ = £(x, t) indicates the displacement of a fluid 
element from the equilibrium position x at the time t, 
whereas the vector Q = V x (£ x B) is the magnetic 
field perturbation. In writing (19), use has been made of 



the fact that the magnetic field curvature, points in the 
direction opposite to that of V/i. 

According to the energy principle (see, e.g., Ref. [351 
[36]), a static MHD equilibrium, is stable if and only 
if SW > for all the allowable displacements £(x, t). 
It is then clear from ( |19[ ), that only two of its terms, 
the second and the fourth, can have a destabilizing ef- 
fect. The second term is associated with kink instabili- 
ties, and represents the interaction of the parallel current 
with the magnetic perturbation. The fourth term, on the 



other hand, accounts for the interchange instability, and 
is destabilizing in the presence of unfavorable curvature, 
that is, when Vp and Vh are antiparallel. A stabilizing 
effect, on the other hand, is due to the work needed to 
bend field lines, which is represented by the first term of 
(19), the energy in the field perpendicular to the equi- 
librium magnetic field. The parallel component of the 
first term of (19), as explained in Ref. [35], also contains 



a stabilizing contribution due to plasma compressibility. 
Finally, a stabilizing contribution due to compressibility 
of sound waves, comes from the third term in (19). 



In order to compare (19) and (18), we must of course 



relate the expression ( 19 1, which is valid for general com- 



pressible ideal MHD, to the specific geometry and order- 
ing of CRMHD. Using the relation Q\\ = Q B/B = 
B(2£ • Vh) + (1/B)$ ■ Vp, obtained from Ref. [35], in the 
limit V • £ = (compressibility effects will be treated 
separately), we obtain 



Q 



B 2 



B 



= (Qj_ + 2 £ • Vh B) 2 

Ql+^-Vh) 2 , (20) 

where Qx = Q Q\\B/B is the perturbation perpendicu- 
lar to the equilibrium magnetic field and, in the last step, 
we assumed that to leading order B = z. The perpendic- 
ular magnetic perturbation, on the other hand, can be 
related to a perturbation of the magnetic flux function 
Sip in the following way: 



Q± = V x (£ x B) = V x (Stpi). 
Thus, expression p0| can be rewritten as 

Q-^b) ~|V^| 2 + 4(C-Vft) 2 



(21) 



(22) 



Neglecting toroidal displacement, so that £-2 = 0, from 
(21 1 one gets also the relation 



6^ z + VT = -£ ■ VV> z + 1 x z, 



(23) 



where T is a scalar function and where, again, we con- 
sidered the toroidal field to be constant at leading order. 
The z component of ( 23 1 then leads to 



Stp 



(24) 



which relates an arbitrary displacement £ to a magnetic 
flux perturbation dip, in the presence of a given equilib- 
rium flux function ip. 

Using the same approximations, the second term of 



( 19 1 can be written as 



- j dxdy x B) • Q ~ - j dxdy j„(£ x B) Q 



= - / <fcdl/.7||(£xB).Qx 



dxdy in (| x B) ■ (VStp x z) 



= / dxdy Sip | • Vjn , 



(25) 
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where, for the last step, we integrated by parts assuming 
V • £ = 0. 

From Eq. (16 1, restricted to the case v(ib) = 0, we 
obtain 



(26) 



which, when inserted into ( 25 1 using d24h , yields 



J dxdy^it xB)-Q 



dxdy 



M 2 ( : F"-p"p/(3-p' 2 /l3-2p"h) 



-2£ • Vp £ • Vft 



(27) 



From (26), on the other hand, one can derive the relation 



T"{iP)-p"^)p(4,)/p-p'\iP)/P-2p"^)h 



x B p 

B 2 



Vj|| 



(z x B p ■ Vp) (zxB p - V/i) 



I? 2 



5 2 

p 



which, when combined with (27) yields 



| dxdy |(£xB)-Q 



dxdy 



Viii 



(28) 



+2 (ixB,V P )(ixB,V fe)) _ 2HpHfe 



_B 2 



B 2 



Finally, by combining Eqs. (19), (22), and (28), assuming 
V • £ = 0, we obtain 



SW 



dxdyi \VSip\ 2 + 4(£- Vhf 



z x B 



B 2 

p 

2 (zxBp-V P ) (z x B p ■ Vfe) 2 



B 2 
p 



I? 2 



£ -Vj|||^| 2 
(29) 



By comparing (29) with (18), we observe that the ex- 



pression for SW contains the same destabilizing terms 
obtained from the EC method. Stability conditions ob- 
tained from the energy principle, are then intimately re- 
lated to those obtained from the general EC method. 
This similarity confirms then that, for the equilibria un- 
der consideration, an indefiniteness in the sign of S 2 F for 
CRMHD, can be attributed to the presence of kink or 
interchange instabilities. 

Of course, CRMHD does not coincide exactly with the 
ideal MHD system to which the energy principle was 
applied, so that an exact one-to-one correspondence be- 
tween SW for ideal MHD and S 2 F for CRMHD should 
not be expected. However, further correspondence can 
be found. First, we observe that the stabilizing field line 
bending term | V<5'0 1 2 appears in both expression. The 
expression for SW contains the additional term propor- 
tional to (£• V/i) 2 , which is however, always positive, and 



therefore, not a possible source for instabilities. Second, 
we remark that for ideal MHD, pressure perturbations 
are related to the displacement by 



Sp 



-i ■ vp - p v • £. 



(30) 



Using ( 24 ) and the fact that p is a flux function, we obtain 



pV-£ = Sp — p 1 (ip)Sip 



= Sp 



x Bp • Vp 
B 2 



Sip. 



(31) 



Consequently, the fourth term under the integral in ( 18 ) 
can be written as 



z x Bp ■ Vp 

B 2 
"p 



Sib 



(32) 



This term is clearly reminiscent of the third term in ( 19 ) 
(recall that (3 in the denominator of S 2 F is proportional 
to the background plasma pressure). This analogy then 
indicates that the stabilizing role of the fourth term in 
S 2 F for CRMHD, is due to the plasma compressibility. 

Finally, we observe that the terms |V<5</>| 2 and \Sv\ 2 in 
S 2 F obviously have no counterpart in SW. Indeed, such 
terms refer to variation of the kinetic energy, and are not 
part of the potential energy SW. 

Now, we extend the EC stability analysis for CRMHD 
to the more general case in which the equilibrium pos- 
sesses also a toroidal flow. Therefore we impose Q = 0, 
but AT ^ 0, so that at equilibrium, v(ip) — —N(ip). 

The second variation, then reads 



5 F 



dxdy 



|V50| 2 + \VS^j\ 2 + \Sv\ 2 + \5p\ 2 /(3 



+ (.F"(V) + {p/P + 2h ) C"(il>) + v Af"(<iP)) |<5</f 



+ 2£'(ib) Sp Sip/ (3 + 2Af'(ip) 8v Sip 



(33) 



Using the equilibrium relations of Eqs. ( 11 )— ( 14 ) with 
Q = 0, the expression (33) can be rewritten as 



5 2 F 



dxdy 



\V8cp\ 2 + |V<5V>| 2 



If z x Bp • Vp 

A 5p m~ H 



z x Bp • Vu 

Sv — Sib 

p 



zxB 



B 2 
"p 



(z x Bp • Vp) (z x Bp • V/i).^, 2 



B 2 

p 



B 2 

p 



(34) 



Comparing Eq. (34) with (18), we observe that, after 



the inclusion of toroidal flow, the only terms with indef- 
inite sign are still those associated with kink and inter- 
change instabilities. Therefore, we conclude that, accord- 
ing to CRMHD, the possible instabilities for equilibria 
with toroidal flow obtained by extremizing the free en- 
ergy functional, can still only be of kink or interchange 
type. 
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III. ANALOGY WITH THE MRI MODEL 



Eqs. (35 )— (38 ) becomes 



Now we consider the model for MRI in the shearing 
sheet approximation derived in Ref. [TBJ. This model 
considers the local behavior of a conducting incompress- 
ible fluid rotating about the z axis of a Cartesian coor- 
dinate system (x, y, z), with x and y indicating the ra- 
dial and azimuthal directions, respectively. The system 
is also supposed to be translationally invariant along the 
azimuthal coordinate, so that all fields depend on x and z 
only. A linear shearflow, Uo = oxy along the azimuthal 
direction, is assumed to be maintained and accounts for 
the radial variation of the angular velocity. The geometry 
of the system is depicted in the left panel of Fig. [T] 

Neglecting dissipative terms, the model equations in a 
dimensionless form, are given by 



dui 

lh 

dip 

~dt 

db 
dt 

do 
m 



+ 2Q[x,v] z + [<j),u}] z -v 2 a [x,W 2 i/j} z 

+ v 2 A [w 2 ip,i/j] z = o, 

+ {<f>,i)] z -[x,<f>] z = Q, 



(35) 
(36) 



+ [<p, b] z + [v, ip] z - [x, v] z + a[x, ip) z = 0, (37) 



+ [(p,v] z + v A [b,ip] z -(2Q + a)[x,i 
-v\ [x,b] z = 0. 



(38) 



In this system ip(x,z) and (p(x, z) denote a magnetic 
flux function and a stream function, respectively, with 
u) = V 2 (p. So, the radial and vertical components of the 
magnetic and velocity fields are given by B x = —dip/dz, 
B z — dip/dx and v x = —dep/dz , v z = dep/dx, re- 
spectively. (Note, in order to facilitate comparison with 
CRMHD, we interchanged the meaning of ip and cp 
adopted in Ref. [IS].) The fields b(x, z) and v(x, z), on 
the other hand, represent the inhomogeneous azimuthal 
components of the magnetic and velocity fields, respec- 
tively. The constant parameters fi, a, and va denote the 
angular velocity, the amplitude of the shearflow and the 
Alfven speed based on a constant background magnetic 
field directed along z. We let [/, g] z = —y ■ V/ x V,g, 
which is the canonical bracket with respect to the vari- 
ables x and z between two functions / and g. 

We transform to the following new variables: 



(p = (p, ip = v A (ip + x), b = v A b, v = v + ax, (39) 



where ip is a flux function that includes a constant back- 
ground z-component of the magnetic field and v is the 
y-component of a velocity field that includes a linear 
shearflow. 



duj 
~di 
dip 
~di 
db 
dt 
di> 
~di 



+ [4>,Q] z + [v 2 ip,ijj} z -2n[ij,x] z -- 

+ [0,6], + [5,^ = 0, 

+ [tp,v] z + [b,i>] z -2n[x,4>] z = o, 



(40) 
(41) 
(42) 
(43) 



where u> = V 2 (p. Comparing Eqs. (40)-(|43j) with the 
CRMHD model of Eqs. Q-Q, it is evident that the 
two models are identical provided one sets /3 = 1 and 
makes the following identifications: 

(p <H- (p, ip -H- ip, b -o- v, v O p, (44) 
y O z, nx^-h, [,]z**[,]y (45) 



It emerges then, from this comparison, that there is an 
analogy between the role played by the angular rotation 
in the MRI system with that of magnetic curvature in 
CRMHD. Indeed, the last terms in Eqs. d40b and d43l, 



which are due to the Coriolis force acting on the fluid, 
match the last terms of Eqs. ([!]) and Q, which come 
from the toroidal component of the Lorentz force (to- 
gether with the assumption that the pressure and the 
diamagnetic correction of the toroidal magnetic field are 
proportional) and the compressibility of the E x B flow in 
the presence of a curved magnetic field, respectively. The 
third term in the pressure equation Q, which also ac- 
counts for compressibility effects, matches the magnetic 
tension term in the azimuthal velocity equation (43). As 
pointed out in Ref. [25], the inclusion of compressibil- 
ity terms (i.e., those proportional to (3) in CRMHD is 
not strictly consistent with the ordering adopted to de- 
rive the model. However, such terms were retained in 
the model due to their "special qualitative importance" 
[25] . It is worth noticing that, without such terms, the 
analogy between the pressure equation of CRMHD and 
the azimuthal velocity equation of the (incompressible) 
MRI model (which does strictly follows an imposed order- 
ing) would not exist. We remark also that the toroidal 
velocity of CRMHD, plays a role analogous to that of 
the azimuthal magnetic field, and in particular, the last 
term of Eq. expressing the parallel pressure gradient, 
mirrors the last term of Eq. ( |42| , which represents the 
parallel gradient of the azimuthal velocity. The analogy 
between the vorticity equation ([!]) and the Ohm's law 
Eq. Q of CRMHD and their counterparts in the MRI 
model is more natural, apart from the already discussed 
curvature / Coriolis terms . 

In spite of this formal analogy between the mod- 
els, however, we remark that, under the transformation 



Eq. ( 39 1 , the new variables do not satisfy the same radial 



In terms of the variables (<p,ip,b,v), the system of boundary conditions as the old. 



IV. 



MRI MODEL HAMILTONIAN STRUCTURE 
AND STABILITY 



For the Hamiltonian structure of the MRI model ( 35 1 



(38), the boundary conditions imposed on the fields are 



analogous to those imposed for CRMHD; that is, vanish- 
ing Dirichlet boundary conditions along x and periodic 
along z. With these assumptions, the Hamiltonian of the 
system is given by the functional 



H = i / dxdz{\\74>\ 2 + (v + ax) 2 



+ v 2 A \Vi>\ 2 + v 2 A b 2 ), 



(46) 



where each term has a clear interpretation: the first two 
terms account for the kinetic energy (including that of the 
shearflow) and last two terms represent the magnetic en- 
ergy. Unlike CRMHD, no internal energy term is present, 
as is typical of incompressible models where the pressure 
is not an independent dynamical variable. 

The Poisson bracket is given by 

{F,G} = J dxdz[u[F u ,G u ] t + b([F u ,G h ], + [F b ,G u ] s ) 

+ (v + (2f> + a)x) ([F u , G v ] z + [F V ,G U ],) 
+ y + x)([F u ,,G i ,] z (47) 

+ [Fjji, G u ] z — [Ft,, G v ] z — [F v , Gt,]z)) , 
which satisfies all the requisite Poisson bracket proper- 



ties. Together with the Hamiltonian (46), this bracket 



yields the MRI model (35)-(38). Moreover, it possesses 
the same structure as the CRMHD bracket of Eq. 
Consequently, the Casimirs of (47), which are given by 



the following four infinite families: 

C\ = J dxdz F{ip + x), 

C 2 = J dxdz b M{tp + x), 

C 3 = J dxdz C(ip + x)(v + (20 + a)x) , 

Ci — J dxdz(uj Q(ip + x) 

-bg'ty + x) (v+(2Ct + a)x)), 

are analogous to those of Eqs. ([8| . We find here, again the 
conservation of magnetic flux, expressed by the family 
Ci, remnants of magnetic helicity, C2, and cross- helicity 
conservation, C3 and C4. 

The EC stability analysis proceeds as for CRMHD. Ex- 
tremizing the free energy functional leads to the following 



equilibrium equations: 

<t> 



v\b 



v + ax — 



G(rf> + x), (48) 
F'(tp + x) + bN''(ip + x) 
+ (« + (2n + a)x)C'{^ + x)+lo g'{tj) + x) 

-bg"(^ + x)((v + (2n + cr)x), (49) 
-7V(-0 + x) 

+ g'(^ + x)( v + (2n + <j)x), (50) 

-£(!/> + x) + b g'ty + x). (51) 



from which a plethora of equilibrium states are possible. 
However, as a first case for this model, we consider the 
trivial equilibrium state 



= 0, ^ = 0, 



b = 0, 



v = 0, 



(52) 



i.e., only the background constant magnetic fields and 
the linear shearflow are present at equilibrium. 

The state ( 52 ) corresponds to the following choice for 



the free functions of the Casimirs (48): 
2VL + a 



T 
C 



-<r(ip + x), 



-(^ + x) 2 , Af=0, 



0. 



The second variation of the corresponding free energy 
functional is then given by 



S 2 F = 



dxdz 



|v<M 2 + H 2 +^iWi 2 

+ v 2 A \Sb\ 2 + (a 2 + 2na)\5iP\ 2 - 2a SipSv 
which can be rearranged as 



S 2 F = 



dxdz 



|V^| 2 + vi|V^| 2 



+ v\\5b\ 2 + \aS1p - 6v\ 2 + 2na\6^ 



(53) 



It is clear from (53) that the only possible nonpositive- 



definite term, that is the only possible source for instabil- 
ity, comes from the last term of (53). Indeed, the second 
variation is not manifestly positive definite if Qa < 0, 
which corresponds (assuming Q > 0) to an angular ve- 
locity radially decreasing outward. This is, in fact, the 
condition for MRI in the shearing sheet approximation 
derived in Ref. [TB]. (Note, application of the Poincare 
inequality for suitable 5^ gives another stability condi- 
tion that depends effectively on the system size.) 

We remark that in CRMHD the trivial state <f> = 
if) = v = p = 0, analogous to (52), corresponds to 
F' = g = JV=C = 0. From (17), one concludes 



immediately that such an equilibrium is always stable. 
Notice, however, whereas the MRI model embodies the 
ingredients for the MRI (identified by the parameters O 
and a, in addition, of course, to the imposed background 
magnetic field) in its equations of motion, this is not the 
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case for CRMHD, which accounts for the field curvature 
but not necessarily for the pressure gradient, which is re- 
quired for the interchange instability. In order to allow 
for this instability to take place, it is therefore neces- 
sary to have a nontrivial solution for the pressure in the 
equilibrium equations (11)-(14). However, apart from 
the trivial states, the analogy between the MRI and the 
interchange instability, already suggested in Sec. be- 



comes evident if one considers more general equilibrium 
states. Indeed, by taking advantage of the analogies be- 
tween the Hamiltonian and the Casimirs in the two mod- 
els, we can readily write the second variation for the free 
energy functional F = H + C\ + C2 + C3 for the case of 
the MRI model. The resulting expression is 



6 F 



dxdz 



\V6<p\ 2 +v 2 A \V5ip\' 



(54) 



Bp 



B 2 



B 2 



2 (y x B p ■ V{v + ax)) (y x B p ■ V{2flx)) 2 



B 2 



B 2 



where B p = V(V> + x) x y and j\\ ■ 
d 2 ip/dz 2 ). It is evident upon comparing (54) with (34 1, 



-{d 2 ip/dx 2 



that the term responsible for the interchange instability 
in CRMHD is analogous to the last term of (54), which 



is responsible for the MRI. If the equilibrium azimuthal 
flow is reduced only to the imposed shearflow (i.e., v = 0), 
then the latter term yields the above mentioned instabil- 
ity condition fla < 0. However, Eq. (54) offers now a 



stability condition that applies to more general equilib- 
rium flows. Also, the next to last term in Eq. (54) indi- 



cates the possibility for a kink-type instability that the 
rotating magnetized disk can undergo in the presence of 
an equilibrium current density. 

Finally, note that if the presence of the magnetic field 
is removed from the MRI model, then the system of 
Eqs. (35H38) reduces the ideal fluid model 



^ + 2Q[x,v] z + [<f>,Lj] z = 0, 
at 

dv 

— + [4>,v] z - 2Q + (T )[x,<f>] s = 0. 
ot 

This model is characterized by the Hamiltonian 



(55) 
(56) 



11 = ^ I dxdz{\V<f)\ 2 + (v + ax) 2 ), (57) 



and the bracket 



{F,G} 



dxdz 



w [F w , G u 



+ (v+ (2fl + a)x) ([F U ,G V ] Z + [F V ,G U ] 9 ) 



which possesses the two families of Casimirs 
Ci = J dxdz F(v + (2fi + a)x) , 
C-2 = I dxdz oj g(v + (20 + a)x) , 



for arbitrary F and Q . The second variation of the free 
energy functional, for this model, evaluated at the trivial 
equilibrium <f> = v = (corresponding to F(v + (2f2 + 
a)x) = -(cr/(4Q + 2a))(v + (2il + a)x) 2 and Q = 0), 
reads 



S 2 F = / dxdz |V<50| 



2tt 



2Q, 



-\5v\ 



(58) 



One can then easily see that, if one assumes, as before, 
51 > 0, energy stability is achieved if 



2n + a > 0, 



(59) 



which is namely the Rayleigh hydrodynamic stability 
condition for the rotating fluid that one recovers from 
Ref. [37]. The regime of relevance for MRI is of course 
that for which the equilibrium is hydrodynamically sta- 
ble (i.e., 251 + a > 0), but the presence of the imposed 
vertical magnetic field makes it energy unstable, which 
occurs when a < 0, as shown above. 



V. CONCLUSIONS 

We have presented analytical results concerning the 
Hamiltonian structure and the stability properties of two 
reduced models, one of interest for tokamak dynamics 
and the other for accretion disks dynamics. Our anal- 
ysis revealed analogies between the two models and, in 
particular, analogous instability processes that they can 
describe. We also demonstrated a connection between 
the standard MHD energy principle and the EC stability 
method. 

More specifically, after reviewing the Hamiltonian 
structure of the CRMHD model, we derived energy sta- 
bility conditions for a wide class of equilibria obtained 
from a variational principle. Two terms of S 2 F that can 
lead to instability were identified. The physical mean- 
ing of these terms was exposed by tracing their mean- 
ing back to the ideal MHD SW analysis, traditionally 
adopted in plasma physics. Thus, these terms were seen 
to corresponding to the kink and interchange instabili- 
ties. Similarly, in the expression for S 2 F, the stabilizing 
roles of compressibility and magnetic field line bending 
were noted. By comparing the cases of equilibria with 
and without toroidal flow, we concluded that instabili- 
ties (such as, e.g., Kelvin-Hclmholtz instabilities) intro- 
duced by the addition of sheared toroidal flow can be 
associated with the terms responsible for the kink and 
interchange mechanisms. Thus, instabilities instigated 
by toroidal flow can then still be detected by considering 
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the direction of current and density gradients, with re- 
spect to magnetic flux gradients and magnetic curvature, 
respectively. 

Next, we described the reduced model of Ref. [TB], in- 
troduced for investigating MRI in the shearing sheet ap- 
proximation. We provided a formal mapping between 
CRMHD and this MRI model. The roles played in the 
latter by the azimuthal magnetic and velocity compo- 
nents, turn out to correspond to those played by toroidal 
velocity and pressure, respectively, in CRMHD. We re- 
marked that, curiously, the analogy between these two 
models is only possible because the compressibility term 
has been kept ad hoc in CRMHD, although this term is 
negligible according to the imposed ordering. 

By taking advantage of the formal analogy between 
the models, we obtained the Hamiltonian structure for 
the MRI model and identified four infinite families of 
Casimir invariants expressing conservation laws. 

We then applied the EC method to the MRI model, 
first for a trivial equilibrium state usually adopted for 
MRI studies. For this equilibrium we recovered the in- 
stability condition Via < 0, derived in Ref. [18 , that 
corresponds to the usual condition of angular velocity de- 
creasing radially outward. We then extended the analysis 
to much more general equilibria. From the analogy 5 2 F 
was immediate; whence we concluded that the destabi- 
lizing (nonpositive-definite) term responsible for MRI is 
analogous to the term responsible for interchange insta- 
bility in CRMHD. Thus, the role played by the magnetic 
curvature in tokamaks is analogous to that played by the 
gradient of the angular velocity in accretion disks. This 



term was seen to give rise to interchange or magnetoro- 
tational instability, respectively, when orientation is in 
an unfavorable direction, with respect to the pressure or 
azimuthal velocity field gradients. The analogy was ex- 
tended further, indicating that accretion disk equilibria 
can undergo instability analogous of the kink instability, 
when the gradient of the parallel current density points 
toward an unfavorable direction with respect to V?/>- 

More generally, the analogy permits the transference 
of knowledge about CRMHD stability results to the ac- 
cretion disk model, and vice versa. Thus the Hamilto- 
nian structure of the models provides a convenient and 
unifying framework for investigating properties such as 
conservation laws and stability conditions for both mod- 
els. 
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